Fidelity susceptibility of one-dimensional models with twisted boundary conditions 



Manisha Thakurathi^, Diptiman Scn^, and Amit Dutta^ 
^ Center for High Energy Physics, Indian Institute of Science, Bangalore 560 012, India 
^Department of Physics, Indian Institute of Technology, Kanpur 208 016, India 

Recently it has been shown that the fidehty of the ground state of a quantum many-body system 
can be used to detect its quantum critical points (QCPs). If g denotes the parameter in the 
Hamiltonian with respect to which the fidelity is computed, we find that for one-dimensional models 
with large but finite size, the fidelity susceptibility xf can detect a QCP provided that the correlation 
length exponent satisfies u <2. We then show that xf can be used to locate a QCP even if u > 2 
if we introduce boundary conditions labeled by a twist angle N6, where is the system size. 
If the QCP lies at g = 0, we find that if A'^ is kept constant, xf has a scaling form given by 
Xp ~ d''^^" fig/d^^") if 6* < 2n/N. We illustrate this both in a tight-binding model of fermions 
with a spatially varying chemical potential with amplitude h and period 2q in which v = q, and in 
a XY spin-1/2 chain in which v = 2. Finally we show that when q is very large, the model has 
two additional QCPs at /i = ±2 which cannot be detected by studying the energy spectrum but are 
clearly detected by xf- The peak value and width of xf seem to scale as non-trivial powers of q at 
these QCPs. We argue that these QCPs mark a transition between extended and localized states 
at the Fermi energy. 

PACS numbers: 64.70.Tg, 03.67.-a, 75.10.Jm 



I. INTRODUCTION 

Quantum phase transitions have been studied exten- 
sively for several years. These are transitions which oc- 
cur in the ground state of a many-body quantum system 
as a parameter g in the Hamiltonian is varied across a 
critical valuei"— ; the ground state typically has differ- 
ent orders on the two sides of the critical point. Several 
measures arising in quantum information theory, such as 
entanglement^"—, entanglement entropy^iii^, Loschmidt 
echoii, decoherencei^ and quantum discordi^"— have 
been used to detect the location of a quantum critical 
point (QCP). A number of reviews have appeared on the 
connections between quantum critical systems and quan- 
tum information theory^—. 

The concept of quantum fidelity {F) has proved to 
be particularly useful for detecting the locations of 
QCPs^^Ti^. We consider a one-dimensional system whose 
Hamiltonian H{g) contains a parameter g such that the 
system is at a QCP when g ~Q. Note that a QCP exists 
only in the thermodynamic limit in which the number of 
sites TV — oo. Consider now a system in which N is finite 
but large 1). Given two ground state wave functions 
|'0o(5 + dg/2)) and IV^olff — dg/2)) at two values of the 
parameter which are separated by a small amount dg (we 
assume that the ground states are non-degenerate), we 
define the fidelity as 

F{g,dg) = \{M9 + dg/2)\M9-dg/2))\. (1) 

First order perturbation theory shows that in the limit 
— > 0, F — 1 — > as {dgY- We therefore define the 
fidelity susceptibility (PS) as 



(Note that our definition of the FS differs by a factor of 
2 from the one given in many other papers^^*^) . The FS 
measures how rapidly the ground state changes with g. It 
turns out that near the QCP at 5 = 0, xf shows a large 
peak even for finite values of N because there are a large 
number of very low-energy states which mix with each 
other in a way which changes rapidly with g. Thus the FS 
is able to detect the ground state singularities associated 
with a quantum phase transition without making explicit 
reference to an order parameter. In the limit g — >■ 0, it is 
found thal^i^-^ 
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where v is the correlation length exponent at the QCP 
(namely, for N ^ 00, the correlation length ^ diverges 
as Igl'" as g 0). Eq. ^ holds only if v < 2; the 
divergence in that expression for small g arises from con- 
tributions from the low-energy (critical) modes. How- 
ever, if v > 2, Eq. ^ is no longer useful for finding 
QCPs because, for small g, the contributions to xf from 
high-energy modes are of the same order as (if v — 2) ot 
dominate over (if ly > 2) the terms of order |.g|'^~^. 

We will show in this paper that even if j/ > 2, we 
can use xf to locate a QCP by using twisted boundary 
conditions labeled by a twist angle N9 and taking the 
limit — > in a particular way. The introduction of 
9 allows us to bring the energy of one particular state 
arbitrarily close to zero; this state then contributes a 
term to xf which scales with 6 as 6""^/^ f{g / 6^^'^). In the 
limit of — >■ 0, a plot of xf versus g clearly shows the 
divergence due to Q~'^^^ at g Q^l" , thereby pinpointing 
the location of the QCP. 

The paper is organized as follows. In Sec. II A, we will 
present a simple argument to show that the presence of a 
twist angle NB, along with a gap proportional to g^ , gives 
rise to a xf which exhibits the above scaling form. In 
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Sec. II B, we indicate how this result may be generahzed 
to higher dimensions. In Sees. Ill A-C, we wih illustrate 
how this works in a tight-binding model of noninteracting 
spinless fermions with a periodic chemical potential with 
amplitude h and period 2g, where q is an integer. This 
model has a QCP at h = where v ^ q; it therefore 
allows us to study the scaling of xp for different values of 
J/. In Sec. Ill D, we will illustrate the scaling of xf with 
6* in a spin-1/2 chain which has a QCP with u = 2. In 
Sec. IV, we will consider the case of large values of q and 
numerically show that xf has a scaling form near h ~ ±2 
with non-trivial power laws. (This study is motivated by 
the observation that the Aubry- Andre model, which has 
a quasiperiodic chemical potential, has transition from 
extended to localized wave functions at h = ±2). In Sec. 
V, we will summarize the results presented in this paper. 



II. FIDELITY SUSCEPTIBILITY FOR A 
TWO-STATE SYSTEM 

In this section we will do a simple calculation to de- 
termine the fidelity of a two-state system. The results 
obtained here will be used in the following sections. 



A. One-dimensional systems 

Consider a two-state problem which is governed by a 
Hamiltonian of the form 



H 



—uk ag'^ 
ag'^ uk 
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This form is motivated as follows. In later sections we 
will consider a one-dimensional model in which ±k will 
denote the momenta of two states close to zero energy 
(measured with respect to the Fermi energy), u will de- 
note the Fermi velocity obtained by linearizing the dis- 
persion near zero energy, g will denote a parameter in a 
many-body Hamiltonian such that the QCP lies at g = 0, 
and a is some constant. The form in Eq. (|1]) will be taken 
to be valid only for values of k much smaller than some 
cut-off, i.e., only for low-energy modes. Since the eigen- 
values of H in Eq. ^ arc equal to ±-\/ u^fc^ + a^g'^'^^ we 
see that the excitation spectrum is given hy E = ±uk at 
the QCP (which would be consistent with the dynamical 
critical exponent of a many-body system being given by 
z = 1), while E = g" for fc = is consistent, for z = 1, 
with the definition of the correlation length exponent v. 

After finding the ground state \ipo{9)) of H in Eq. (U), 
we can compute the FS using Eqs. (dill); we find 



form 

XF{9,k) = {uk/a)-'^-- f{a'/''g/{uky^n, 

1/2 1 

where f(x) = ^ ■ — -77. (6) 

We note that f{x) has a single peak at x = if j/ = 1. 
li ly > 1, f{x) has two peaks at a; = ±[(i^ — l)/(j^ + 
and goes to as -{i^^ /8)x^'-''-^'^ for x ^ and 
as -t/V(8x2(''+i)) for x ^ 00. 

Next, let us suppose that we are considering a many- 
body system with size N, so that the momenta are quan- 
tized as k = 9 + 2T:n/N, where n runs over all inte- 
gers, and 9 comes from a twisted boundary condition 
to be elaborated in later sections. The presence of the 
27Tn/N term implies that we can take 9 to lie in the range 
[0,27r/A^]. If the Hamiltonian of the system can be de- 
composed into independent two-state systems of the form 
given in Eq. the FS will be given by 
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We can now consider the value of xf in various limits. 

First, if N is held fixed and we consider a regime in 
which ag" , u9 < 2t:u/N and a}/" g / (119^^" is of order 1, 
the sum in Eq. ([7]) will be dominated by the n = term 
and we will obtain the scaling form 



Xf{9) = {u9/a)-'-/^ f{a}/^g/{u9f/'^), 



(8) 



where the function / is given in Eq. Secondly, if 

^ 2ttu/N but 9 is of order 2tt/N , then we obtain 



Xf{9) 



n— — 00 



8u2 

32u2 s\n^{N9/2)' 



[9 + 2TTn/Nf 



(9) 



As 5 — ^ 0, Eq. (O shows that xf has a finite limit \iv = \ 
but it goes to if > 1. [For the models that we will 
consider later, Eq. (|4]) is not a good approximation for 
|fc| g'^ (high-energy modes), and the contributions to 
the FS from such modes generally approaches a non-zero 
value as (7 —s- as we will show in Eq. for v > \\. 
Finally, if 27r/A^ ^ g", the sum in Eq. ([7]) can be replaced 
by an integral over k = 21:111 N] if we further assume that 
(7" ^ 1, the limits of the integral can be taken to be ±00 
since the contributions from the regions with \k\ ^ g'^ 
will be negligible. We then obtain 



XF{g,k) = - 



1 



(5) 



8g2 [ag''/{uk) + uk/i^ag")]^ 
In the limit fc 0, we see that the FS has the scaling 



Xf{9) = - 



89' .-00 
Ni^^ag"^^ 
32u 



dk 



27r [ag''/{uk) 
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ufc/(ag'')]2 
(10) 
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[The integral in the above equation only gets substan- 
tial contributions from values of |fc| lying around g'^ , i.e., 
from low-energy modes. The scaling given in Eq. (fTO|) 
will therefore be valid as long as Eq. is a good ap- 
proximation for such low-energy modes, even if it fails 
for high-energy modes]. We thus see from Eq. (fTO]) that 
as 5 — > 0, the FS will diverge ]iv<2 but not liv >2. A 
plot of the FS versus g will therefore show a divergence 
at the QCP {g = 0) only if j/ < 2, but for j/ > 2, the FS 
will not be useful for finding the location of the QCP. 

We conclude that for a large but finite value of the 
system size TV, the FS can detect a QCP only if < 2. 
Note that this conclusion would not change if we replaced 
k and g'^ by fc^ and g'^'^ respectively in Eq. Q, with a 
dynamical critical exponent z which is not equal to 1. 
The exponent which appears in the second line in Eq. (|10p 
will be — 2 regardless of the value of z. 



g^ ^ 2ttu/N but is of order 2tt/N, then we obtain 



B. Higher dimensional systems 

Although the rest of this paper will be only concerned 
with one-dimensional systems, let us briefly discuss what 
may happen in higher dimensions. We assume that there 
is a d-dimensional system in which the modes with mo- 
menta ±fc are governed by a Hamiltonian of the form 
given in Eq. (|4]), where the momentum variable in that 
equation now stands for k = \k\. Eqs. ([SKS]) will continue 
to hold. Next, let us impose periodic boundary condi- 
tions, with a twist angle in one of the d dimensions and 
zero twist angle in the remaining d — 1 dimensions. To 
be specific, we assume that we have a hypercubic lattice 
system which has N sites in each dimension, so that the 
different momenta are quantized as ki = 2TTni/N; the 
twist shifts the first momentum to fci ^ 6 + 2TTni/N. 
Eq. ([7]) then gets modified to 



XF 



2 



rii ,rid = -oo 



2 ' 



ag'^ I uk' \ 
uk' ag^ J 

d 

k' = [(0 + 2^ni/iV)2 + ^(2^7i,/Ar)2]i/2. (11) 



As before, we can now consider what happens in three 
different cases. 

If N is held fixed, ag" , uO < 2ttu/N and g/iuOY/" 
is of order 1, the sum in Eq. ([TT]) will be dominated by 
the term in which = for all i. Then we will obtain 
the scaling form given in Eq. ([S]) where the function / 
is given in Eq. Note that the scaling form does not 
depend on the dimensionality d in this case. Next, if 



Xf{9) 
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[NO + 2T:n^Y + YH=2i'^^^i) 



(12) 



As g 0, Eq. (fT2|) has a finite limit (which depends on 
NO) liv^l but goes to ff > 1. Finally, ff 27r/A^ < g" , 
the sum in Eq. (jlip can be replaced by an integral over 
ki = 2TTni/N; assuming that g'^ <^ \, the limits of the 
integral can be taken to be ±oo. We then obtain 



XF{g) 



- N'^a'^g'^"-^ 



d'^k 



1 



W [ag'^/iu\k\) + u\k\/iag'^)]^ 

(13) 



(A similar scaling relation was found by adiabatic per- 
turbation theory^^i^^) . As 5 — >■ 0, Eq. ([T^ diverges if 
dv <2 but not if dv >2. A plot of xf therefore diverges 
at the QCP if dv < 2, but for dv > 2, xf is not useful 
for locating the QCP. 



III. TIGHT-BINDING MODEL WITH 
VARIABLE 

We will now study fidelity in a one-dimensional tight- 
binding model of spinless fermions in which the exponent 
V can take any integer value; this will illustrate many of 
the points discussed in Sec. II. In Sec. Ill A, we will 
discuss the model and its energy spectrum close to zero 
energy, while in Sec. Ill B, we will examine a number of 
features related to fidelity. 



A. The model 

The model that we are interested in was studied re- 
cently from the point of view of quenching dynamics 
across a QCP— 1^. In this section, we will summarize 
the relevant discussion from Ref. [SSl The Hamiltonian 
is given by 

N 
n=l 

+ h C0S(^ + (j)) clcn], (14) 

where g is a positive integer, N is the system size, and 
we have imposed periodic boundary conditions so that 
cn+i = ci. (We will set the hopping amplitude J, h 
and the lattice spacing a equal to unity). Note that one 
can perform a unitary transformation on the c„ , namely, 
c„ — )■ e~'"^c„, which removes the phase 9 from all the 
hopping terms except for the hopping between the sites 
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at A'' and 1 where the phase becomes e^'^-^^; this is called 
a twisted boundary condition, with NO being the twist 
angle. We can assume that < N9 < 2tt, namely, that 
lies in the range [0,27r/A^]. The usual periodic bound- 
ary condition corresponds to 6* = 0. The period of the 
chemical potential in Eq. (fT4|) is 2q^ and we will assume 
that is a multiple of 2q. A chemical potential of this 
form appears in the Azbel-Hofstadter model^"— . Since 
shifting n — > ri + 1 in Eq. (|14p is equivalent to shifting 
(/)—>■ + Tr/g, we can assume without loss of generality 
that (j) lies in the range [0, 7r/g]. 

The fermionic operators can be Fourier transformed to 
the momentum basis. 



Ck 



N 

E 



-ikn 



(15) 



where the momentum k goes from — tt to tt in units of 
2-K /N] these operators satisfy the anti-commutation rules 
{cfc, c\i} = Sk.k' ■ III momentum space, the first two terms 
of the Hamiltonian in Eq. ([T4)) have the tight-binding 
form 



Hn = - 



^ 2cos(fc + 6) 4cfe. 



k——7r 



(16) 



For the last term in Eq. ([T4| . we use the decomposition 

hcos(— + (P) = (e^(W9+0) 4. g-*(W9-i-0)). (17) 
q 2 

Hence this term couples two fermionic modes with mo- 
menta ki and ^2 if fci = ^2 ± Tr/g. This fragments the 
total Hamiltonian H into N/ {2q) decoupled Hamiltoni- 
ans Hk which are labeled by a momentum k lying in the 
range — tt to — 7r + 7r/g, namely, 



-7r+7r/g 

H = Hk. 

k—~1T 



(18) 



For each fc, Hk can be written as a (2(j)-dimensional 
matrix involving the momenta k + rn/q^ where r = 
1, 2, • • • , 2q. The matrix elements of Hk are given by 



(fc + riT/q\Hk\k + sir/q) 
= —2 cos(fc + 6 + sn/q)Sr,! 
-{h/2){e'^Sr,s+i+e-'^Sr, 



(19) 



where 1 < r, s < 2q, and we have set J = 1. In writing 
Eq. (|19p. we have assumed 'periodic boundary conditions' 
for the matrix Hk , so that r ~ and 2q+l mean r = 2q 
and 1 respectively. 

Before proceeding further, let us point out an inter- 
esting consequence of the decoupling of H into a sum 
of Hk and then a duality symmetry^i^. Suppose that 
the system has a parameter 9 and A'^ = 2qr sites, where 
r is a positive integer. Since the momentum k is quan- 
tized in units of 2-k/N , and only momenta differing by 



integer multiples of 7r/g are coupled to each other, the 
form of Hk in Eq. (|19p implies that this system is exactly 
equivalent to a sum of r different systems, which have 
A^' = 2q sites each but have r different values of 9' given 
by 61, 6'+27r/(2gr), 9+AT:/{2qr), ••• , 6'-h2(r-l)7r/(2gr). 
Therefore, instead of studying a system with 2qr sites and 
one particular value of 9, we can study a system with only 
2q sites but several values of 9. Next, we observe that if 
A^ ~ 2q, the Hamiltonians in real and momentum space, 
Eqs. and ([TO)) , get mapped into each other if we si- 
multaneously interchange J ^ h/2 and 9 ^ (p; here we 
have used the fact that a shift in or by Tr/g = 2tt/N 
has no effect on the spectrum. In particular, the model 
in Eq. is self-dual if J /i/2 and 6* 0. 

For each of the Hamiltonians Hk in Eq. ([T^ . the 2q 
energy levels come in q pairs ±E. This can be shown by 
shifting s — > s + g which flips the sign of the first term in 
Eq. ([T5)) , and changing the sign of the state corresponding 
to s by (—1)''' which flips the sign of the terms propor- 
tional to h in that equation. Assuming that there are no 
states at exactly zero energy, we see that the ground state 
of each Hk is one in which the q negative energy states 
are filled and the q positive energy states are empty; thus 
the ground state is half-filled for each Hk- 

It will soon become clear that the model defined above 
has a QCP at /i = 0. We therefore use perturbation the- 
ory to study the region near h = 0^. (We will set 9 = 
in Eqs. (|20II23P : this is a reasonable approximation if A^ 
is large since we always take < 9 < 2tt/N). We will be 
particularly interested in states near zero energy which 
will contribute the most to the FS. Since the system is at 
half- filling, these states lie near the momenta k = ±/cf, 
where the Fermi momentum kp = Tr/2. If the amplitude 
h of the chemical potential is zero, the system is gapless 
and the states aX k = ±7r/2 are degenerate with each 
other. Assuming that h is small compared to the band 
width of 4J, we can use a perturbative expansion in h 
to calculate the breaking of this degeneracy. The regions 
near k = — 7r/2 and k = tt/2 regions are coupled through 
one series of intermediate states lying at fc = — 7r/2-|-7r/g, 
— 7r/2 -|- 27r/(7, • • • , 7r/2 — Tr/g (with an amplitude equal to 
— he^"^ at each step), and through another series of inter- 
mediate states lying at fc = — 7r/2 — -k / q, — 7r/2 — 27r/q, 
• • • , 7r/2 -I- ir/q (with an amplitude equal to — /le^*"^ at 
each step). Each of these series consists of q — 1 inter- 
mediate states. At the g-th order in perturbation theory, 
we therefore obtain an effective Hamiltonian Hg^f f which 
has a matrix element between the states at k = ±7r/2 
given by 



AE^(V2|i/e//|-7r/2) 
{-h/2Y e"J't' 



-V2|i?e//k/2) 



Ill=l {2cos{-n/2 + s7T/q)) 

{-h/2Y e-"i'^ 

^ UlZl (2cos(-V2-W'Z))' 

where the denominators come from factors like E^^ 
E_^/2±sTr/q = 2cos(— 7r/2-|-±s7r/g) corresponding to the 



(20) 
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energies in the unperturbed Hamiltonian in Eq. (1161 
Simplifying Eq. (|20p gives 



{-ir{h/2)'> 



r,iqct> 



UlZl (2sin(Wg)) 
Hence the magnitude of A is given by 



(21) 



hi 


|cos(g0)| 




_J sin(7rs/g) 


hi 


|sin(g0)| 


4,-1 Yii: 


sin{ns/q) 



|A| = -r-— r ' if q is odd, 

if q is even. (22) 



We see that (f) governs the relative phase between the 
two sets of intermediate states which connect the states 
at A: = ±7r/2. We will assume that (j> is such that 
cos((j'(/)) 7^ if q is odd, and sin(g0) 7^ if <; is even; 
if these conditions are violated, we would have to go to 
higher order perturbation theory to find a non-zero ma- 
trix element connecting the states at = ±7r/2. 

We can now consider moving slightly away from k = 
±7r/2; then the unperturbed energies of the states k = 
— 7r/2 + k' and TT/2 + k' are given by —2k' and 2k' respec- 
tively. The effective Hamiltonian describing these two 
states is then given by the 2x2 matrix 



eff,k' 



-2k' A* 
A 2k' 



(23) 



where we assume that A continues to be given by the ex- 
pression in Eq. (PT|) because k' is small. The eigenvalues 
of are given by ± ^4^'^ + |Ap; this is the dispersion 
of a massive relativistic particle whose velocity is equal 
to the Fermi velocity u = 2 and mass is proportional to 
|A| ~ h'^ times cos(g(/)) or sin((70). 

Hence, h = corresponds to a QCP where the mass 
gap vanishes. Given that the energy vanishes as \k'\ if 
h = and as h"^ if k' ~ 0, the dynamical critical exponent 
and correlation length exponent are given by z = 1 and 
v — q, respectively. The correlation length exponent v 
thus depends in a simple way on the periodicity of the 
chemical potential. 

Comparing Eq. (g]) with Eqs. (|22ll23p . we identify g = 
h, u ^ 2, and the constant a is given by 



1 



\cos{q(f>)\ 



4''"' ml sm{ns/q) 

1 \sm{q(j))\ 
4-^-1 UlZl Sinins/q) 



if q is odd, 

if q is even. (24) 



To summarize, we began with a model whose Hamil- 
tonian in momentum space consists of a sum of {2q)- 
dimensional Hamiltonians. Close to the QCP which lies 
at h = 0, we used perturbation theory to write an effec- 
tive two-state Hamiltonian which governs pairs of low- 
energy states, i.e., states close to the Fermi energy. In 
the next section, we will use the results obtained in Sec. 
H A for the FS of two-state systems to obtain the FS of 
our model. 



B. Fidelity susceptibility for various v 

We begin this section by describing how the FS can be 
numerically calculated for the model presented in Sec. 
HI A. As we have seen, the Hamiltonian decouples into 
N/{2q) Hamiltonians Hk- If we can compute the fidelity 
Fk{h,dh) = liipo^kih + dh)\ipo^k{h))\ for each of the Hk, 
the total fidelity of the system will be given by the prod- 
uct 



-TT+TT /q 



F{h,dh) 



Fk{h,dh). 



(25) 



k — — Tl 



Next, we know that the ground state is half- filled for 
each of the Hk for every value of h. Let ipi.k{h) de- 
note the first quantized wave functions of the filled states 
i {i = 1,2, ■ ■ ■ ,q); each of the ipi^kih) denotes a (2*7)- 
dimensional column which is the space on which Hk 
acts. Let 'ipi^r,k{h) denote the r-th component oi ipi^kih) 
(r = 1, 2, • • • , 2q), and ^ denote the second quantized 
creation operator for the fermion corresponding to the r- 
th component. Then the i-th filled state can be written 
in second quantized form as 'ipl j.{h)\vac) , where 



2, 



(26) 



and \vac) denotes the vacuum state of the fcrmions: 
Cr.k\vac) = for all r and k. The half-filled ground state 
is therefore given by the second quantized expression 



^i.ki^)\vac) 



(27) 



i=l 



We can now use Wick's theorem^ to show that 
Fk{h,dh) = |(wac|(g) '4),,k{h + dh/2) 

1 



^ i,\^{h-dhl2)\vac)\ (28) 
j=i 

is given by the magnitude of the determinant of a q- 
dimensional matrix Mk{h, dh), 

Fk{h,dh) = \det[Mk{h,dh)]\, 
{Mk{h,dh))ij = {vac\'ijji^k{h + dh/2)'ijj],,{h-dh/2)\vac) 
2, 

= E r^,r,k{^ + dh/2)i^,.r,k{h-dh/2). 



r=l 



(29) 



Thus the computation of the total fidelity F reduces to 
finding the determinants of the matrices Mk and then 
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multiplying the determinants over N/{2q) values of k. 
The FS is then found as 



XFih) 



lim 

dh-¥Q 



\nF{h,dh) 



(30) 



We can now use the results obtained in Sec. II to 
understand various properties of the FS of the model 
discussed in Sec. Ill A; we have seen that v = q for 
this model. Further, the parameters in Eq. ^ are given 
hy g = h, u = 2, and a is given in Eq. ([St]). If 27r/iV < 
h"^ ^ I, the FS should indicate the location of the QCP 
a q — I, but not if g > 2. 




FIG. 1: (Color online) Plot of X-F versus h for q = 1, (j> ~ 
(red solid), q ~ 2, <j} = n/A (black dashed), g = 3, (p = (blue 
dotted), and q = 4, (j> — n/S (magenta plus), with A'^ = 240 
and e = Tv/N. 

Fig. [1] confirms the above statement for g = 1, 2, 3 
and 4, iV = 240 and 6 = tt/N. For each value of q in 
that figure, we have chosen (j) so as to maximize the gap 
|A| which is proportional to cos{q(j)) (sin(g0)) for q odd 
(even). At /i = 0, there is a large peak in ~xf for 9=1, 
but not for g = 2, 3 and 4. For 9=1, the peak value of 
—Xf is found to be 450 at /i = which lies far outside the 
range of Fig.[Tl we note that this peak value is consistent 
with Eq. ^ since v = 1, a = 1, and the Fermi velocity 
is u = 2. 

In producing Fig. [TJ we have chosen A'^ to be a multiple 
of 4 and 6 = tt/N for the following reason. As discussed 
in Eq. (|23|) . the momentum k introduced in Sec. II is 
actually the deviation from the momenta ±7r/2 for the 
model discussed in Sec. Ill A. To prevent the expression 
in Eq. ([9|) from diverging in our model, wc must ensure 
that ±7r/2 + 6* + 2TmlN does not vanish for any value 
of n. This will be true if A is a multiple of 4 and 9 = 
tt/N. (Alternatively, we could have chosen A^ — 2 to be 
a multiple of 4 and 9 = 0). 

We observe in Fig. [1] that xf is not zero at = 
for q = 2, 3 and 4, in contrast to the result in Eq. 
We can analytically compute xpih = 0) by using the 



following result from first order perturbation theory. If 
H{X) = Hq + XV is a many-body Hamiltonian with eigen- 
values Ea{X) and eigenstates ^paiX), where a = denotes 
the ground state, then the FS is given by^ 



Xf(A) 



2 ^ 

Q#0 



\{i^^{x)\v\Mm' 

{E^{X) - Eo{X)y 



(31) 



We now apply this result to our model. At /i = 
and 9 = tt/N, the ground state of the total system 
is one in which the A^/2 one-particles states with k = 
-7r/2, -7r/2 + 2tt/N, ■■■ ,tt/2- 2tt/N arc filled and the 
remaining A^/2 states are empty. Next, at h = dh, the 
change in the ground state to first order in dh is given by 
a state in which a fermion has moved from a filled state 
at k to an empty state at fc — Tr/g if fc lies in the range 
[— 7r/2, —tt/2 + TT/q], or from fc to fc + Tr/g if fc lies in the 
range [tt/2 — TT/q,TT/2]. Using Eqs. ()29II30|) . we find that 
the FS at /i = is given by 



XF{h = 0) 

-Tr/2+7r/<j 

- E 

fc=-7r/2 



1 



16 [cos(fc + 9)- cos(fc - TT/q + 9)Y' 



(32) 



where k goes in steps of 2tt/N, and we have used 
the fact that the contribution to xf from the range 
[—tt/2, —tt/2 + TT/q] is equal to the contribution from 
the range [7r/2 — TT/q, tt/2]. For large A^, we can change 
the summation in Eq. ([32]) to an integral {J2k ~^ 
J dk{N/2TT)) and ignore 9 to obtain 



XFih = 0) 



N 



32tt sm{TT/q) 



(33) 



Note that this result does not depend on the phase (j) in 
Eq. HH). For A^ = 240 and g = 2, 3 and 4, this gives 
XF{h = 0) ~ 2.39, 2.76 and 3.38 respectively, which agree 
well with the values shown in Fig. [T] 

Finally we note in Fig. [T] that the FS for q = 2, 3 and 
4 has a double peak structure; however, these peaks are 
unrelated to the QCP at /i = 0. We notice that as q 
increases, these peaks move towards h = ±2. We will see 
in Sec. IV that in the limit g — > oo, there are QCPs at 
h = ±2 in addition to the QCP at h — 0, and the double 
peaks are an indication of those QCPs. 



C. Scaling of fidelity susceptibility with 6 

In this section, we will show that by varying 9, we can 
locate the QCP at /i = even for g > 2. Further, the 
scaling of the FS with respect to 9 allows us to find the 
value of the critical exponent ly. At the end we will also 
discuss the scaling of the FS with (j>. 

In Fig. H wc show the FS for g = 2, 3 and 4, iV = 240 
and 9 = 0.001, which is much smaller than the value of 
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h 



FIG. 2: (Color online) Plot of xf versus h ior q = 2, (j} = tt/A 
(black dashed), q = 3, (j> = (blue dotted), and g = 4, <f> — 
Tv/8 (magenta plus), with N = 240 and 6 = 0.001. 



9 = tt/N chosen in Fig. [TJ the values of 4> chosen for each 
q are the same as in Fig. [TJ We now see that additional 
double peaks have appeared in the different curves which 
lie much closer to h = 0, with the peak values decreasing 
and their distances from h = increasing as q increases. 
We will now see that these new peaks are related to the 
QCP at /i = and that their peak values and positions 
scale in accordance with Eq. ([5]); as 6* — >■ 0, the locations 
of the peaks approach the QCP. 



240 and = 0, with 6 = 0.0001, 0.00003 and 0.00001. 
Once again the curves fall on top of each other. 

The values of 9 in Figs. [3] and |4] indicate that as q 
increases from 2 to 3, we have to go to smaller values of 
9 to see the scaling form. The reason for this is as follows. 
We saw below Eq. ^ that the peak in the FS occurs at 
g/{u9y^i = [{q- l)/((7 + l)]i/(2«) ^hen c/«,6' < 2nu/N. 
In our model, g = h, u = 2, and a is given in Eq. (|24p . 
Putting these together, we find that the peak occurs at 
9 = (a/3/8) h"^ for g = 2 and a.t 9 = (a/2/24) for q ^ 3. 
We thus see that for a small value of h, 9 is smaller for 
g = 3 compared to q = 2. For h = 0.1, for instance, the 
peak lies at 6* ~ 0.0022 for g = 2 and at 6* ~ 0.000059 for 
q = 3. 
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FIG. 3: (Color online) Plot of xfO versus /i/f^/^ for 9 = 
0.0004 (black dashed), 9 = 0.0002 (blue dotted), and 9 = 
0.0001 (magenta plus), with g = 2, iV = 240 and (j) = 7r/4. 

In Fig. [21 we show xf^*^^' versus h/9^/'i for q = 2, 
N = 240 and </> = 7r/4, with 9 = 0.0004, 0.0002 and 
0.0001. We see that the three curves fall on top of each 
other, thus confirming the scaling form given in Eq. ([5|). 
In Fig. [1 we show xpd^^'' versus h/9^/'' for g = 3, iV = 



FIG. 4: (Color online) Plot of xf9^^^ versus h/9^^^ for 9 = 
0.0001 (black dashed), 9 = 0.00003 (blue dotted), and 9 = 
0.00001 (magenta plus), with g = 3, iV = 240 and <?!> = 0. 



— 



2 



FIG. 5: (Color onhne) Plot of xf9 versus h/9^^'^ for = 
(black dashed), (fi — tt/S (blue dotted), and 4> = tt/A (magenta 
plus), with g = 2, = 240 and 9 = 0.0004. 
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FIG. 6: (Color online) Plot of xfO versus /i| sin(2<;/>)|^/V^^''^ 
for = 7r/16 (black dashed), = 7r/8 (blue dotted), and 
(j) = 7r/4 (magenta plus), with g = 2, iV = 240 and = 0.0001. 



It is interesting to consider the dependence of the FS 
on the phase </>. At /i = 0, we saw that xf is given by 
Eq. ([55)1 and is independent of 0. But xf does depend 
on (j) away from /i = 0. Fig. [S] shows xfS versus 
ioT q = 2, N = 240 and 6 = 0.0004, with = 0, 7r/8 
and 7r/4. The three curves look quite different; in fact, 
the curve for = hardly changes with within 
the range shown in the figure. Note that </> = is the 
value at which the perturbatively calculated gap A given 
by Eq. (|22p vanishes for g = 2 and is therefore inde- 
pendent of h within the pcrturbative range. Hence xf 
is essentially independent of h and is simply given by 
its value at ft, = if </> = 0. Using Eqs. (O and (|24l) . 
we find that ior q = ly = 2, xfO should be a function 
of h\sm{2(j))\^/^/9^/^. In Fig. [HI we show xf^ versus 
h\ sin(20)|i/V6'i/2 for q = 2, iV = 240 and 9 = 0.0001, 
with (f) = 7r/16, 7r/8 and 7r/4. We see that the three 
curves fall on top of each other thus confirming the pres- 
ence of the factor of | sin(20)|-'^/^ in the scaling function. 



D. Fidelity susceptibility in another model with 

u = 2 



H 



E 



{Jx + Jy) (XX , „v „y 



i^J-x Jy) 



(ft-(-l)"<5) 



(34) 



The spectrum of this model can be solved by carrying 
out the Jordan- Wigner transformation from spin-l/2's to 
spinless fermions at each site. (The system then decou- 
ples into a number of four-dimensional subsystems. We 
will not present the details here and refer the reader to 
Refs. m and [13). Defining + Jy — J and Jx — Jy = 7, 
we find that the model has four quantum critical lines 
given by S ■ 



± a/Zi^ -I- 7^ and h 



±V5^TJ^. If 5 = ±7 
is held fixed, there are QCPs at /i = which lie on the 
critical lines S = ±y7?+7^. Alternatively, if ft. — ±J 
is held fixed, there are QCPs at (5 = which lie on the 
critical lines ft = \fW^+lP' . All these QCPs have z = \ 
and V ^2. 

We have numerically studied the FS of this model af- 
ter introducing a twist Q in the Jordan- Wigner fermionic 
Hamiltonian. Holding J = 2, 7 = 1 and h ~ 2 fixed, we 
varied 5 to go through the QCP lying at 5 = 0. Fig. [7] 
shows the scaling of xf^ versus S/e^/^ for N ^ 240 and 
9 = 0.0004, 0.0002 and 0.0001. The scaling near the QCP 
at (5 = is very similar to what is seen in Fig. |3l just as 
one would expect from Eq. (|5]) for a system with v = 2. 
This model therefore confirms that the introduction of 
a twist angle can enable the FS to detect a QCP with 
1^ = 2. 




It turns out that there is another model in one dimen- 
sion which has QCPs with v = 2. This is an anisotropic 
spin- 1/2 XY spin chain in which the strength of a trans- 
verse field alternates between h + S and ft — (5 at odd and 
even sites^"— . The Hamiltonian of the model is given 



FIG. 7: (Color online) Plot of Xff versus 5/0^/^ for 6 = 
0.0004 (black dashed), 9 = 0.0002 (blue dotted), and 6 = 
0.0001 (magenta plus) for the anisotropic spin-1/2 chain with 
N = 240. 
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IV. FIDELITY SUSCEPTIBILITY FOR VERY 
LONG PERIODS 

If the quantity ir / q in Eq. (|14p is replaced by tt times 
an irrational number (which can be approximated by ra- 
tional numbers with increasingly large denominators), we 
obtain a quasiperiodic systen*^"— i^i^. It is known that 
such a system has a metal-insulator transition at ft, = ±2; 
the nature of the eigenstate at zero energy changes from 
extended (metallic) to localized (insulating) on crossing 
these QCPsi^i^i^. The FS has been used to detect 
QCPs in quasiperiodic systems^ as well as in disordered 
systems^. 

In this section, we will study the FS of the model de- 
fined in Eq. ([H]) as q becomes very large. In that limit, 
we find numerically that the FS has increasingly large 
peaks at /i = ±2; we have seen precursors of these peaks 
in Fig. 1 for q = 2, 3 and 4. As far as we know, these 
QCPs have not been reported earlier. Our results in- 
dicate that the metal-insulator transition also seems to 
occur in our model, although 1/q approaches zero rather 
than an irrational number as g — > (X). 
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FIG. 8: (Color online) Plot of xf versus hiox q = 120 (black 
dashed), q — 240 (blue dotted), and q — 480 (magenta solid), 
with N = 2q, 6 — tt/N and (p = n/N in each case. 

In Fig. [51 we show the FS for q = 120, 240 and 480, 
with N = 2q, 9 ^ tt/N and (f> = n/N in each case. (Since 
N must be a multiple of 2q and q is quite large, we are 
only presenting the results for N = 2q here. However, we 
have checked for g = 60 that our results do not change if 
we take iV to be a higher multiple of 2q. Further, we have 
fixed the values of 6 and 4> in terms of A^, and will not 
study here how the FS varies with them). We observe 
prominent peaks in the FS at /i = ±2. Fig. ^ shows a 
log-log plot of the peak value and the full width at half 
maximum versus q for the peak in xf at h = 2. for a 
range of values of q from 30 to 2280. We find that the 
peak value of xf scales as g2.25 -^j^jig ^^j^g f^u width at 
half maximum scales as l/g"'^". In Fig. [TOl we show the 
scaling near ft, = 2 by plotting Xf/*?^'^^ versus (ft— 2)(7° ™ 
for q ~ 120, 240 and 480; the curves fall on top of each 
other. Interestingly, we see that the curves are somewhat 
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FIG. 9: (Color online) Log-log plot of the peak value (blue 
solid, y-axis on left) and the full width at half maximum (red 
dashed, y-axis on right) versus q for the peak in xf at ft = 2. 
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FIG. 10: (Color online) Plot of xf/(7^-^^ versus (/i-2)g" ™ for 
q = 120 (black dashed), q = 240 (blue dotted), and q = 480 
(magenta solid), with N = 2q, 9 = tt/N and = -n/N in each 
case. 



asymmetric about ft = 2. 

We have also studied the energy gap and the nature of 
the wave function at zero energy, i.e., the Fermi energy. 
Upon extrapolating to the thermodynamic limit N ^ oo, 
we find that the gap is zero even if we are away from 
ft ^ ±2; hence these QCPs are different from the ones 
discussed in the earlier sections where the gap scales as 
I ft — hc\'^ and is therefore non-zero for h ^ he even if 
N ^ oo (here he denotes the location of the QCP). The 
QCPs at ft = ±2 are not characterized by the gap going 
to zero but rather by a change in the nature of the wave 
function at zero energy. 

We can understand why the behavior of the system 
changes at ft = ±2 by using a continuum theory to study 
the properties of Eq. for a state whose energy lies at 
zero, i.e., at the Fermi energy. (We will assume here that 
the limit N ^ oo has been taken). A continuum theory 
is justified if q is very large since the chemical potential 
then varies on a length scale which is much longer than 
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the lattice spacing a. We can remove the twist angle 9, 
by performing the phase transformation c„ — > e~'"^c„ in 
Eq. (fT4|)). if we are only interested in the behavior of a 
state in a local region of the system. Setting J = 1, the 
equation of motion following from Eq. (jl4p is 

- (c„+i + c„_i) - /icos( h 0) c„ c„ (35) 

for a state with energy i?. If c„ varies slowly with n, 
we can write 2c„ — c„+i — c„_i = '~cPc{x) / dx^ , where 
X = na (we will set a = 1). Assuming /i > 0, we rede- 
fine n as n — TiQ where no is an integer chosen in such 
a way that ttuo /q + (f> is a.s close to tt as possible; hence 
— hcos(Tmo/q + (f>) is close to its maximum value of h. 
(We can choose uq in an infinite number of ways; the 
various choices differ from each other by multiples of 2q) . 
We then find that Eq. ([55]) can be written as a differential 
equation for c{x), 

d"^ C ITT 

- ^4 + [/icos( — ) -2\c = Ec. (36) 
ax'^ q 

Eq. p6p describes a particle moving in a periodic poten- 
tial whose maximum value is h — 2. The potential has 
an infinite number of wells, each extending from uq to 
no -t- 2q. If /i < 2, a state with E ^ lies above the max- 
imum of the potential; in that case the particle can move 
classically between the different wells, and the wave func- 
tion will be extended throughout the system. If /i > 2, a 
particle with i? = is classically confined to one particu- 
lar well of the periodic potential and can only go to other 
wells by quantum mechanical tunneling. If h is slightly 
greater than 2, a WKB approximation shows that the 
tunneling probability between two neighboring wells is 
proportional to exp[— (ft, — 2)q]. If {h — 2)q >> 1, the 
tunneling probability is extremely small, and the wave 
function is localized within a single well. We thus see 
that in the limit (7 — >■ oo, ft = 2 marks a transition be- 
tween extended and localized wave functions for the state 
with zero energy. 

V. CONCLUSIONS 

To summarize, we have shown that in a one- 
dimensional model which has a QCP with a correlation 
length exponent v > 2, the introduction of a twist angle 
9 enables us to use the fidelity susceptibility to determine 
the location of the QCP. Namely, if the Hamiltonian has 
a parameter g such that the QCP lies at g = 0, the twist 
allows us to bring the energy of a particular state close to 
zero. If 61 < 2n/N, the FS scales as xf - d'^^" fig/d^^") 
which makes the QCP clearly visible if xf is plotted ver- 
sus g. A twisted boundary condition therefore provides 
a powerful tool for locating a QCP which may be diffi- 
cult to find in any other way. We have argued that this 
technique may also be useful for finding QCPs in higher 
dimensional models. 



The specific model that we have used to demonstrate 
this idea is a tight-binding model of spinless fermions 
with a periodic chemical potential with amplitude ft and 
period 2q, where q is an integer. We have studied in 
detail the QCP lying at ft = 0; at this point, the dy- 
namical critical exponent is given by z = 1 while v = q. 
This makes this model specially useful for testing what 
happens for different values of v. We have shown an- 
alytically, using perturbation theory about ft = and 
the decoupling of the system into a number of {2q)- 
dimensional subsystems, that the fidelity susceptibility 
scales as xf ~ 9~^/'^ f(h/9^^'^) if 9 is sufficiently small; 
we have verified this scaling numerically for small values 
of q. Although we have not presented the details here, 
we have confirmed that a similar scaling relation holds 
in another model which has QCPs with v — 2. This is 
an anisotropic spin- 1/2 XY chain with transverse fields 
which has alternating strengths on odd and even sites. 

In the last part of the paper, we have considered what 
happens in our model when q becomes very large. We 
find that some additional QCPs appear at ft = ±2 in that 
limit, and the FS is clearly able to detect these QCPs. 
To the best of our knowledge, these QCPs have not been 
reported before. We have studied the power laws associ- 
ated with the peak value and the width of the FS as a 
function of g; we find non-trivial powers of 2.25 for the 
peak value and 0.70 for the width. The energy gap be- 
tween the ground state and the first excited state is zero 
in the thermodynamic limit N 00 over a finite range 
of values of ft around these QCPs. This makes it diffi- 
cult to define the exponent v at these QCPs (unlike the 
QCP lying at ft = where we know that v = q). Using a 
continuum theory, we have argued that these QCPs are 
characterized by a change in the nature of the wave func- 
tion of a particle at the Fermi energy from extended to 
localized. In the future it would be interesting to develop 
a more detailed understanding of these QCPs which may 
shed some light on the non-trivial power laws that we 
have found. 

The analysis in this paper is expected to be valid for 
any one-dimensional system which reduces to a theory 
of non- interacting fermions in the low-energy limit, i.e., 
close to the Fermi energy. Near a QCP with z = 1 and an 
arbitrary value of the correlation length exponent v, the 
modes near the critical momenta (= ±7r/2 for a half-filled 
system) can be described by 2 x 2 Hamiltonians as in Eqs. 
(HI and p3| . From this we can deduce the scaling of the 
fidelity susceptibility with g and 9. We also argued at 
the end of Sec. II A that the analysis can be generalized 
to the case of z not equal to 1 and in Sec. II B to theories 
of non-interacting fermions in higher dimensions. 

The case of interacting systems is more complicated. 
In one dimension, such systems are typically described 
by Tomonaga-Luttinger liquid theory which has z = 1. 
In a recent paper—, some of us studied quenching in an 
interacting system which is closely related to the model 
considered in this paper. We argued there that the ef- 
fect of interactions is to change the value of v from q 
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to q/{2 — K) where the parameter K depends on the 
strength of the interactions, provided that K Ues in the 
range 1/2 < K < 2. (A non-interacting system has 
if = 1, and we then recover the resuh v = q). We expect 
that a similar resuh would hold for the fidelity suscepti- 
bility of interacting systems; this may be an interesting 
question for future investigation. 

Finally, we note that the twisted boundary condition 
has been introduced in this paper as a mathematical de- 
vice to produce a non-trivial scaling of the fidelity suscep- 
tibility which can provide information about the value of 
v. However, such a boundary condition also admits an 
interesting physical interpretation. A one-dimensional 
system with periodic boundary conditions is the same as 
a circle. Imposing twisted boundary conditions in such 
a system is equivalent to assigning a charge to the par- 



ticles and passing a magnetic flux through the middle of 
the circle so that the Aharonov-Bohm phase (the product 
of the charge and the magnetic flux in some appropriate 
units) is equal to the twist angle N9. By introducing a 
twist wc arc therefore effectively studying the effect of a 
magnetic flux on the fidelity susceptibility of a system of 
charged particles moving on a circle. 



Acknowledgments 

For financial support, M.T. and A.D. thank CSIR, 
India and D.S. thanks DST, India for Project No. 
SR/S2/JCB-44/2010. 



^ S. Sachdev, Quantum Phase Transitions (Cambridge Uni- 
versity Press, Cambridge, 1999). 

^ B. K. Chakrabarti, A. Dutta, and P. Sen, Quantum Ising 
Phases and Transitions in Transverse Ising Models, Lec- 
ture Notes in Physics, Vol. m41 (Springer, Heidelberg, 
1996). 

^ M. A. Continentino, Quantum Scaling in Many-Body Sys- 
tems (World Scientific, Singapore, 2001). 

* S. L. Sondhi, S. M. Girvin, J. P. Carini, and D. Shahar, 
Rev. Mod. Phys. 69, 315 (1997). 

^ M. Vojta, Rep. Prog. Phys. 66, 2069 (2003). 
® A. Osterloh, L. Amico, G. Falci, and R. Fazio, Nature 
(London) 416, 608 (2002). 

T. J. Osborne and M. A. Nielsen, Phys. Rev. A 66, 032110 
(2002). 

* L. Campos Venuti, C. Degli Esposti Boschi, and M. 
Roncaglia, Phys. Rev. Lett. 96, 247206 (2006). 

® G. Vidal, J. I. Latorre, E. Rico, and A. Kitaev, Phys. Rev. 

Lett. 90, 227902 (2003). 
^° A. Kitaev and J. Preskill, Phys. Rev. Lett. 96, 110404 
(2006). 

" H. T. Quan, Z. Song, X. F. Liu, P. Zanardi, and C. P. Sun, 
Phys. Rev. Lett. 96, 140604 (2006). 

B. Damski, H. T. Quan, and W. H. Zurek, Phys. Rev. A 
83, 062104 (2011). 
" H. OUivier and W. H. Zurek, Phys. Rev. Lett. 88 017901 
(2001). 

" R. Dillenschneider, Phys. Rev. B 78, 224413 (2008). 

S. Luo, Phys. Rev. A 77, 042303 (2008). 

M. S. Sarandy, Phys. Rev. A 80, 022108 (2009). 
" T. Nag, A. Patra, and A. Dutta, J. Stat. Mech. (2011) 

P08026. 

L. Amico, R. Fazio, A. Osterloh, and V. Vedral, Rev. Mod. 
Phys. 80, 517 (2008). 

J. I. Latorre and A. Rierra, J. Phys. A 42, 504002 (2009). 
A. Polkovnikov, K. Sengupta, A. Silva, and M. Vengalat- 
tore. Rev. Mod. Phys. 83, 863 (2010). 
A. Dutta, U. Divakaran, D. Sen, B. K. Chakrabarti, T. F. 
Rosenbaum, and G. Aeppli. „arXiv: 1 012.0653^2 (2010). 
M. Znidaric and T. Prosen, J. Phys. A 36, 2463 (2003). 
P. Zanardi and N. Paunkovic, Phys. Rev. E 74, 031123 
(2006). 



L. Campos Venuti and P. Zanardi, Phys. Rev. Lett. 99, 
095701 (2007). 

P. Zanardi, P. Giorda, and M. Cozzini, Phys. Rev. Lett. 
99, 100603 (2007). 

P. Buonsante and A. Vezzani, Phys. Rev. Lett. 98, 110601 

(2007) . 

W.-L. You, Y.-W. Li, and S.-J. Gu, Phys. Rev. E 76, 
022101 (2007); S. Yang, S.-J. Gu, C.-P. Sun, and H.-Q. Lin, 
Phys. Rev. A 78, 012304 (2008); S.-J. Gu, H.-M. Kwok, 
W.-Q. Ning, and H.-Q. Lin, Phys. Rev. B 77, 245109 

(2008) ; S.-J. Gu and H.-Q. Lin, EPL 87, 10003 (2009); 
S.-J. Gu, Int. J. Mod. Phys B 24, 4371 (2010). 

H.-Q. Zhou, R. Orus, and G. Vidal, Phys. Rev. Lett. 100, 
080601 (2008); H.-Q. Zhou, J. H. Zhao, and B. Li, J. Phys. 
A 41, 492002 (2008); H.-Q. Zhou and J. P. Barjaktarevic, 
J. Phys. A, 41 412001 (2008); J.-H. Zhao and H.-Q. Zhou, 
Phys. Rev. B 80, 014403 (2009). 
^3 J. Ma, L. Xu, H.-N. Xiong, and X. Wang, Phys. Rev. E 
78, 051126 (2008). 

E. Eriksson and H. Johannesson, Phys. Rev. A 79, 
060301 (R) (2009). 

D. Schwandt, F. Alet, and S. Capponi, Phys. Rev. Lett. 
103, 170501 (2009); A. F. Albuquerque, F. Alet, C. Sire, 
and S. Capponi, Phys. Rev. B 81, 064418 (2010). 
^2 V. Gritsev and A. Polkovnikov, arXiv:0910.3692 (2009), 
published in Understanding Quantum Phase Transitions, 
edited by L. D. Carr (Taylor and Francis, Boca Raton, 
2010). 

C. De Grandi, V. Gritsev, and A. Polkovnikov, Phys. Rev. 
B 81, 012303 (2010); C. De Grandi, V. Gritsev, and A. 
Polkovnikov, Phys. Rev. B 81, 224301 (2010). 

M. M. Rams and B. Damski, Phys. Rev. Lett. 106, 055701 
(2011); M. M. Rams and B. Damski, Phys. Rev. A 84 
032324 (2011). 

J. Sirker, Phys. Rev. Lett. 105, 117203 (2010). 

M. Thesberg and E. S. S0rensen, Phys. Rev. B 84, 224435 

(2011). 

^'^ V. Mukherjee, A. Polkovnikov, and A. Dutta, Phys. Rev. 
B 83, 075118 (2011); V. Mukherjee and A. Dutta, Phys. 
Rev. B 83, 214302 (2011); V. Mukherjee, A. Dutta, and 

D. Sen, Phys. Rev. B 85, 024301 (2012). 

M. Thakurathi, W. DeGottardi, D. Sen, and S. Vishvesh- 



12 



wara, Phys. Rev. B 85, 165425 (2012). 
^® W. DeGottardi, D. Sen, and S. Vishvesliwara, New J. Phys 
13, (2011) 065208; D. Sen and S. Vishveshwara, EPL 91, 
66009 (2010). 

*° M. Ya. Azbel, Zh. Eksp. Teor. Fiz. 46, 929 (1964) [Sov. 
Phys. JETP 19, 634 (1964)]; D. R. Hofstadter, Phys. Rev. 
B 14, 2239 (1976). 

W. Y. Hsu and L. M. Fahcov, Phys. Rev. B 13, 1595 
(1976). 

^■^ S. Aubry and G. Andre, Ann. Israel Phys. Soc. 3, 133 
(1980). 

G. M. Obermair and H.-J. Schellnhuber, Phys. Rev. B 23, 

5185 (1981); H.-J. Schellnhuber, G. M. Obermair, and A. 

Rauh, Phys. Rev. B 23, 5191 (1981). 

S. Ostlund and R. Pandit, Phys. Rev. B 29, 1394 (1984). 

J. B. Sokoloff, Phys. Rep. 126, 189 (1985). 

F. Delyon, J. Phys. A 20, L21 (1987). 

S. N. Sun and J. P. Ralston, Phys. Rev. B 44, 13603 (1991). 

P. B. Wiegmann and A. V. Zabrodin, Phys. Rev. Lett. 72, 

1890 (1994). 



D. Sen and S. Lai, Phys. Rev. B 61, 9001 (2000), and 
Europhys. Lett. 52, 337 (2000). 

C. Aulbach, A. Wobst, G.-L. Ingold, P. Hanggi, and I. 
Varga, New J. Phys. 6, 70 (2004). 
M. Modugno, New J. Phys. 11, 033023 (2009). 
C. Itzykson and J.-B. Zuber, Quantum Field Theory 
(McGraw-Hill, Singapore, 1980). 
'^^ S. Deng, G. Ortiz, and L. Viola, EPL 84, 67008 (2008). 
U. Divakaran, A. Dutta, and D. Sen, Phys. Rev. B 78, 
144301 (2008). 

^■"^ J. H. H. Perk, H. W. Capel, and M. J. Zuilhof, Physica 

81A, 319 (1975). 
"'^ K. Okamoto and K. Yasumura, J. Phys. Soc. Jpn. 59, 993 

(1990). 

" S. Wen-Guang and T. Pei-Qing, Chin. Phys. B 18, 4707 
(2009). 

''^ J. C. C. Cestari, A. Foerster, and M. A. Gusmao, Phys. 
Rev. A 82, 063634 (2010). 



